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We make an estimate of the likelihood function for the Higgs vacuum expectation value by im- 
posing anthropic constraints on the existence of atoms while allowing the other parameters of the 
Standard Model to also be variable. We argue that the most important extra ingredients are the 
Yukawa couplings, and for the intrinsic distribution of Yukawa couplings we use the scale invari- 
ant distribution which is favored phenomenologically. The result is successful phenomenologically, 
favoring values close to the observed vev. We also discuss modifications that can change these con- 
clusions. Our work supports the hypothesis that the anthropic constraints could be the origin of 
the small value of the Higgs vev. 
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INTRODUCTION 



It is known that if the masses of the hght quarks 
and the electron were modestly different then nuclei and 
atoms would not exist [ll-0]- Because the masses of 
fermions are proportional to the Higgs vacuum expec- 
tation value (vev), these bounds can be interpreted as 
constraints on possible values of the Higgs vev if the other 
parameters of the Standard Model are held fixed [l|, Q . 
This observation is interesting because it could provide 
an answer to one of the most significant puzzles of the 
Standard Model - often called the fine tuning problem 
or the hierarchy problem. If the underlying theory al- 
lows the existence of different values of the Higgs vev, we 
would only find ourself in a region of the universe that 
contains atoms, and hence the vev may be constrained 
to a small range of possible values. This provides a mo- 
tivation for theories with multiple possible values of the 
physical parameters, such as the string theory landscape 
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However, in theories in which the basic parameters 
can take on multiple values, other parameters besides 
the Higgs vev will most likely also be variable. This 
would be the case in the string landscape picture. Since 
the atomic constraints are really on the up quark, down 
quark and electron masses, they translate to constraints 
on the product of the fermion Yukawa couplings and the 
Higgs vev, not the vev uniquely. Even without the explo- 
ration of specific theories, we might hope that the rough 
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conclusion is unchanged, namely that the scale of the 
electroweak sector must be reasonably close to the scale 
of the strong interactions in order that the masses of the 
light quarks and electron - the product of the electroweak 
interactions - be comparable to nuclear binding energies 
- primarily due to strong interactions. 

In this paper we consider the more general case of al- 
lowing other parameters of the Standard Model to vary, 
and attempt to provide a likelihood distribution for the 
Higgs vev. One might think that this would be possible 
only with the knowledge of the full underlying theory, 
but we will primarily use data for this purpose. As we 
will argue in Sec. 2 and 3, this is possible because it is the 
Yukawa couplings that appear to have the most signifi- 
cant infiuence on the range of the vev. Because there are 
many masses, we then have experimental indication of 
the intrinsic probability distribution for the Yukawa cou- 
plings. The observed quark and lepton masses provide 
quite strong statistical evidence that this distribution is 
close to scale invariant Q. We will review this idea 
in Section 4 and provide further evidence in its favor in 
Section 5. 

Applying such a scale invariant probability distribution 
for the Yukawa couplings, we investigate the likelihood 
distribution for the Higgs vev. This result is obtained by 
finding the relative probability that the u, d quarks and 
the electron (governed by the scale invariant weight) fall 
in the anthropically allowed range. In this case our result 
is developed and displayed in Sec. 6. We can also study 
the effect of producing modest changes in our underlying 
assumptions. These are studied in Sec 7. We present our 
conclusions in Sec 8. 

We are aware of many limitations of our work. Be- 
sides the assumptions that we state and explore, there are 
likely other effects (nucleosynthesis, cosmology, etc) that 
come into play, especially once we consider significant 
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changes in the parameters of the Standard ModeL How- 
ever, one would expect that possible further anthropic 
constraints would only tighten the likelihood function in 
the neighborhood of the physical value. The goal of the 
present work is to obtain a sense of whether or not the 
atomic constraints could be the origin of the low value of 
the Higgs vev (within the context of landscape-like theo- 
ries). Our work can be viewed as an attempt to quantify 
this by looking at what may be the dominant effects. 
Overall our conclusion is that it remains plausible that 
the atomic constraints are the origin of the low value of 
the Higgs vev. 



2. GENERAL FRAMEWORK 

The situation that we have in mind is similar to the 
string landscape picture in which there are very many 
possible values of each of the parameters. While in string 
theory the choices of parameters are discrete, the results 
appear to be so densely packed as to appear almost con- 
tinuous^ . We then describe the ensemble of such states by 
an intrinsic probability distribution or weight that spec- 
ifies the probability of finding different values of the pa- 
rameters. These probabilities would emerge from string 
theory and the weight encodes the shape of the string 
landscape. Let us call this weight p{v,Ti,gi) where v is 
the Higgs vacuum expectation value, Ti are the Yukawa 
couplings, and gi stands for the gauge couplings and all 
other parameters of the theory. 

However, many combinations of the parameters do not 
lead to nuclei and atoms. There is an intrinsic selection 
effect that we would only find ourselves in friendly regions 
that include atoms. The shape of the intrinsic probabil- 
ity distribution in unfriendly regions of parameter space 
is then completely irrelevant for us and we are only con- 
cerned with the parameter subspace that leads to atoms. 
Let us denote by A{v,Ti, gi) the function that is zero for 
all parameters that do not lead to atoms and unity for 
those that do. We will refer to this as the atomic function. 
In principle, the atomic function could also take into ac- 
count not only the mere existence of atoms but also the 
probability of a physical environment of sufficient com- 
plexity developing with the atoms that are available with 
that parameter set. As one moves around the parame- 
ter space, especially near the allowed borders, greater 
or fewer numbers of atoms exist and/or would be pro- 
duced in the early universe. With a reduced or enhanced 
set of atoms available, the resulting complexity might be 
greatly reduced or enhanced. However, such considera- 
tions are beyond our capabilities to calculate and we do 



For example, it has been estimated that there are 10^°° string 
vacua reproducing the Standard Model parameters within the 
present experimental error bars Q, and the density of states 
would be equally high in the neighborhood of these parameters. 



not consider them. Given the primary features uncov- 
ered below, it is unlikely that modifying the boundaries 
of the atomic function would have a large effect on our 
results. Moreover, it is certainly cleaner and more con- 
servative to limit our discussion to the general physical 
characteristics of atoms and nuclei. 

With the atomic function we can obtain the total prob- 
ability to find atoms in the landscape 

P{A) = J d.v dTi dgi A{v,ri,gi) p{v,Ti,gi) (1) 

where A denotes the existance of atoms. However, this 
total probability is not a quantity of interest. A more 
interesting one is constructed by omitting the integration 
over V, 

L{v) = j dVi dgi A{v,T.i,gi) p{v,r,,gi) (2) 

which we call the likelihood function. It is in fact the 
probability density ^ for atoms to exist. 

Another useful assumption is the independence of pa- 
rameters. This means that we assume that the intrinsic 
probability function factorizes into the product of sepa- 
rate weights. In formulas this implies 

p{v, Lj, g,) = p{v)p{re)p{Tu)p{Td) (3) 

This is at least partially motivated by the vastness of 
the string landscape. If we hold all but one parameter 
fixed, there are likely other allowed vacua with this last 
parameter scanning over its allowed range. There is also 
some phenomenological evidence in favor of this from the 
distribution of quark and lepton masses, which all seem 
consistent with the same distribution. However, we note 
the approximate nature of this feature in our discussion 
of Sec. 7. 

There also could be an a-priori distribution for the 
Higgs vev, which is a property of the fundamental 

theory. We clearly do not know this function. However, 
it is expected that the vev can take on values in a very 
large range, at least up to a unification scale. Moreover, 
there are many additive contributions to v that come 
from quantum corrections. These add linearly in the 
respective couplings and this suggests that the overall 
distribution could be Taylor expanded in v about the ob- 
served value. If this is the case, then our considerations 
cover only a very small portion of the allowed range, and 
we treat the a-priori distribution as a constant in this 
narrow range. If the a-priori distribution in v were to 
be highly peaked in some direction, our results would be 
modified, and so this must count as an uncertainty in our 
method. Therefore, our assumption will be p{v) ^ con- 
stant, and unless we know more about the underlying 
theory, we feel that this is the most reasonable assump- 
tion under which to proceed. 

The quantity we will explore in detail below is the 
probability to find atoms for a given value of It is 
obtained by taking a sample v and drawing the other 
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relevant parameters randomly from the probability dis- 
tributions we consider. Then we decide if the resulting 
configuration can yield atoms or not. With a large sam- 
ple for a fixed value of v we can obtain the probability 
of having atoms by dividing the number of times we ob- 
tained atoms by the total number of simulations. Using 
the assumption of independence of the parameters intro- 
duced above, the quantity we obtain from our simulations 
then is 

P(A|given v)=J dT, J dg, [x{p{Vi)p{gj)^ A{T,,g,,v). 

(4) 

Now the likelihood function L{v) of Eq. ([2]) is obtained 
by the product of P(A|given v) and the intrinsic proba- 
bility distribution p{v). Under our assumption of a fiat 
p{v) ~ const, in the range of interest, the likelihood func- 
tion is then simply proportional to P(A|given v) of Eq. 

©• 

From the shape of the likelihood function L{v) we can 
infer which values of v are typical and which ones are 
highly improbable. Since L{v) is a probability density 
its shape itself is not a direct indicator of the most likely 
values of v. A peak in L{v) for example does not indicate 
the most likely values of w; more meaningful quantities to 
give would be the median or other percentiles. A simpler 
way to explore the order of magnitude of the most likely 
values of v can be obtained by plotting our results for 
L{v) on a log-log scale. Since any probability distribution 
which has a finite value of its percentiles must fall of 
faster than 1/v aX large values of w, the log- log plots show 
us if and when the likelihood function falls off faster than 
1/v. If present, this point is then a reasonable estimate of 
the most likely values of v. If L{v) does not fall off faster 
than l/u, no constraints on the Higgs vev arise from the 
existance of atoms. 



3. BRIEF SUMMARY OF ATOMIC 
CONSTRAINTS 

To the extent that we understand how the Standard 
Model leads to the world that we observe, we should be 
able to describe the world that would result if we instead 
used parameters different from, but in the neighborhood 
of, those seen in Nature. Surprisingly, the structure of the 
elements changes dramatically for quite modest changes 
in the quark masses. In a recent paper Damour and 
Donoghue have tightened and summarized the anthropic 
constraints on quark masses^ . Here we briefly summarize 
these results. 

The first constraint which results from the binding of 
nuclei gives an upper bound on the sum m„ -f md. The 
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FIG. 1: The anthropic constraints on m^, rud, me in MeV 
units. 



key feature here is that the pion mass-squared is pro- 
portional to this sum of masses, and as the pion mass 
gets larger nuclear binding quickly becomes weaker. The 
binding energy is small on the scale of QCD and is known 
to have opposing effects from an intermediate range at- 
traction and a shorter range repulsion. The attractive 
component, heavily due to two pion exchange, is the most 
sensitive to the pion mass and weakening it leads to a lack 
of binding of nuclei. From [i| this constraint is 

mu + md< 18 MeV. (5) 

The second constraint comes from the stability of pro- 
tons. If protons could annihilate with electrons, p+e^ — > 
n + hydrogen would not exist. The proton and neu- 
tron mass difference gets contributions from the quark 
masses and from electromagnetic interactions. Using the 
best present estimates of these, the constraint becomes 

Oi] 

md-mu~ 1.67me > 0.83 MeV. (6) 

The right hand side of the equation is linear in the elec- 
tromagnetic fine structure constant. Modest variations 
in this number would not infiuence our results signifi- 
cantly. In providing this constraint, it has been assumed 
that the neutrino masses remain negligibly small. This 
feature is also anthropically required J7j . 

These constraints are summarized in Fig. [TJ Note that 
the up quark and electron masses are able to vary down 
to zero mass, while the down quark mass is constrained 
to be non-zero. The dimensional scale is set by the QCD 
scale Age d , so that these constraints could be rephrased 
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in terms of dimensionless ratios rrii/ Kqcd^ ■ 

There are no known atomic constraints on the masses 
of the heavier quarks and leptons as long as they are sig- 
nificantly heavier than the up, down and electron. Heavy 
quarks decouple from low energy physics and have little 
influence on nuclei and atoms. Therefore, in our case the 
atomic function A{v^Ti,gi) reduces to A(t;, r„, r^, Fe). 
However, the Higgs vev does influence the mass of the W 
gauge boson, and this in turn influences the rates of weak 
processes. Since weak interactions play a significant role 
in the pp cycle in stars, in particular the crucial reaction 
p-\-p — )■ d + e^ -\-v, raising the W mass will slow the pro- 
duction of elements. So while atoms may exist, they may 
not be produced in favorable numbers. This constraint 
is complicated also by the dependence on the cosmologi- 
cal parameters governing stars. In 8], a special scenario 
was constructed in which the W mass can be taken to 
infinity if modifications are made to the usual picture of 
nucleosynthesis^. So while there may be additional con- 
straints from the W mass which would further tighten 
the likelihood function, they have not been convincingly 
shown to be tighter than those of the atomic function, 
so we will ignore them in our analysis. Because further 
constraints lead to a narrowing of the likelihood function, 
omitting them is therefore a more conservative approach. 



4. THE DISTRIBUTION OF QUARK AND 
LEPTON MASSES 

In a landscape picture, the Yukawa couplings would 
not be uniquely determined but would follow from some 
intrinsic distribution. Because there are enough masses, 
the observed masses can give us insight into what this 
distribution is without having to know the full underlying 
theory of the landscape [1, [g^- 

The intrinsic probability distribution for a quark or 
lepton Yukawa coupling is defined such that the fraction 
of values that appear at coupling T within a range dV is 
p{T) dV , with the normalization 
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j dV p{T). 



(7) 



In particular, we have explored a set of power law weights 
ii 

(8) 



p(^^ I rmin <^ r 



otherwise, 



^ In fact, A.QCD serves as the comparison scale for all dimensional 
quantities in this work, so that the reference value for the Higgs 
vev Vf) is also in units of A.qcd ■ 

* See, however, Q which raises problems with this scenario. 

^ The Yukawa interactions also influence quark mixing and the 
observed weight is consistent with the hierarchy of weak mixing 
elements |qjl0|[- There are also possible implications for neutrino 
properties [q. Illl . 
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FIG. 2: Quark and lepton masses, defined at the energy 
fi = Mw, on a logarithmic scale. A scale invariant weight 
corresponds to a uniform distribution on this scale. We use 
this feature to describe the distribution of Yukawa couplings 



where the normalization constant is = (1 — ^)/[rJjjai — 
r]-,^] if <5 ^ 1 and iV = 1/ ln[r„ax/r„in] if 5 - 1. 

Such simple power law weights require at least one end- 
point (or two for (5 = I) in order to be normalizable as 
in Eq. ([?])■ When we determine the endpoints for the 
Yukawa distributions at large and low values, it is natu- 
ral to use the renormalization group quasi-fixed point [12| 
^max 1-26 as the upper limit. In [1, 0] a lower end- 
point Tmin ~ 1.18 X 10~^ which corresponds to OAnie 
was used. This is explored more in the following section. 
With these ingredients, a likelihood analysis found that 
(5 = 1.02 ±0.08. 

Of particular interest is the experimentally favored 
scale invariant distribution corresponds to a weight with 
,5 = 1: 



(9) 



While the evidence for the scale invariant weight is based 
on quantitative studies, the result can be seen qualita- 
tively in Fig. [2l A scale invariant distribution is one 
which is a random uniform population on a logarithmic 
scale. Fig. [5] shows the Yukawa couplings for the quarks 
and leptons, at the scale p = Mw, plotted on a logarith- 
mic scale. The result appears visually to be consistent 
with this idea, and in practice a scale invariant weight is 
highly favored^. We will use this as our primary weight 
for our analysis. 

In exploring the uncertainties in the Higgs likelihood 
function, we will also consider weights which have no 
lower bound on the Yukawa couplings. This is only pos- 
sible for 5 < 1, if the distribution is to be normalizable. 



® Because the result is consistent with scale invariant, the choice 
of scale is not important and a similar result would be obtained 
using either the un-rescaled masses or with the Yukawa couplings 
defined at the Grand Unification scale. 
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FIG. 3: The staircase-like curve shows the empirical distri- 
bution function, defined as the fraction of the 9 observed 
Yukawa couplings that lie below a given F-value. The other 
curves show the cumulative distribution functions predicted 
by our power-law probability distribution Ansatz. All assume 
Tmax = 1-2, and the straight lines correspond to the log- 
uniform 5=1 case with different values of Fmin, the solid 
line having the best-fit value Fmin = 5 x 10~^. From top to 
bottom, the bent curves with the same endpoints take the 
5-values 1.2, 1.1, 1.0, 0.9 and 0.8, respectively. The 6 = 1 line 
is seen to exchibit good consistency with the observed data. 

While the possibility Tmin = is statistically disfavored 
(see next section), we found in ^] that the best fit value 
in that case is 5 = 0.86tao5- We will use this in our tests 
of uncertainties in Sec 7. 



5. FURTHER EVIDENCE ON THE FERMION 
MASS DISTRIBUTION 

Since the Bayesian result of Rcf. S = -1.02 ± 0.08 
from a likelihood analysis does not address the ques- 
tion of whether the best-fit distribution is actually con- 
sistent with the data, we have studied a set of fre- 
quentist Kolmogorov-Smirnov (KS) tests. The obser- 
vational input into a KS-test is the so-called empiri- 
cal distribution function J^obs(r), defined as the frac- 
tion of the 9 observed Yukawa couplings that lie be- 
low a given F-value. This function is plotted in Fig. [3l 
and the horizontal location of the staircase-shaped steps 
correspond to the 9 rank-ordered Yukawa couplings for 
(e, /i, T, u, d, s, c. h, t), taken from the Particle Data 
Group compilation [l3| . For any particular choice of our 
model parameters [5, Fmin, Tmax), we can compute a pre- 
dicted cumulative probability distribution function 

p(F')dF' (10) 



FIG. 4: The probability that our model is consistent with 9 
observed Yukawa couplings is shown as a function of the power 
law index 5, assuming the best-fit value Fmin = 5 x lO"'^. The 
V-shaped curve shows the maximal difference between the 
predicted and observed distribution functions from Fig. [S] 



and compare how well it agrees with i^obs(r)- Fig- El 
illustrates this for variations in both 5 and Fmin- The 
KS-test uses as a goodness-of-fit statistic the maximum 
vertical discrepancy [Fobs(r) — ^'(r)| between theory and 
observation. The corresponding consistency probababil- 
ity is shown in Figs. |4] and [5] for variations in i5 and Fmin, 
respectively, keeping a fixed upper cutoff Fmax = 1-2. 

The constraints on 5 from Fig. |4] nicely reproduce 
the previous likelihood result using frequentist methods. 
Fig. |4] shows that we can reject the null hypothesis that 
our model is correct at high significance if 5 deviates 
substantially from unity, and also that the data are per- 
fectly consistent with 5 = 1. This approach also en- 
ables us to place a bound on the lower endpoint Fmin, 
as seen in Fig. [SJ Normalizability alone requires a lower 
endpoint when 5 > 1, but for any 5- value, it is clear 
from Fig. [3] that the fit becomes poor if the left end- 
point is dragged sufficiently far to the left of the left- 
most data point. Fig. [5] shows that the best fit value is 
Fmin ~ 5 X 10~^ for the scale-invariant case, although the 
constraints are rather weak, with 1 — cr and 2 — a bounds 
(32% and 0.045% consistency probability) corresponding 
to Fmin « 7 X 10~^ and w 5 x 10~^^, respectively. 

In Ref. a Bayesian approach was employed where 
the dependence of the likelihood function on the expo- 
nent 5 was studied for a fixed Tmin ■ Here we extend this 
analysis by investigating also the dependence of the like- 
lihood on Train ■ As in [6| we fix the upper endpoint to be 
Fmax = 1.26 which is motivated by the quasi- fixed point 
of the Standard Model. In Fig. |6] we show the contour 
plot of the log-likelihood function as a function of the 
lower Yukawa endpoint Train and the exponent 5. The 
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FIG. 5: The probability that our model is consistent with 9 
observed Yukawa couplings is shown as a function of the lower 
cutoff Fmin, assumiug the scale-invariant power law index S — 
I. The V-shaped curve shows the maximal difference between 
the predicted and observed distribution functions from Fig. [S] 



darkest area is the 1 — a range, the next one is the 2 — a 
range etc. Here the n — a range is taken as the parameter 
space where the log-likelhood is at most below its 

maximum. 

We see that for all values of 6 the data favors a lower 
endpoint for the distribution, and for any fixed value of 5 
the likelihood function increases monotonically up to the 
highest possible value Tmin = Te- The best fit, i.e. the 
highest likelihood, is found for i5 = 1.06 and Tmin = Te- 
Any probability distribution with Tmin = Te seems of 
course very unnatural since one out of the nine measured 
Yukawas would lie exactly on the endpoint of the prob- 
ability distribution, but the likelihood analysis does not 
take that into account. For S < 1 where the power law 
weights do not require a lower endpoint Fmm, the data 
says that such a scenario with Fmm = is disfavored by 
over 2 ~ a. 



6. THE LIKELIHOOD FUNCTION FOR THE 
HIGGS VEV 

We now combine our two key ingredients, the an- 
thropic constraints and the probability distribution for 
the Yukawa couplings which is phenomenologically suc- 
cessful in describing the observed Yukawa couplings, in 
order to estimate the likelihood distribution for the Higgs 
vev. Our approximation of the full problem, under the 
assumptions described in Section 2, consists of 



0.7 0.8 0.9 1.0 1.1 1.2 1.3 1.4 

5 

FIG. 6: Contour plot of the log-likelihood as a function of the 
lower endpoint for the Yukawa couplings Tmin and the ex- 
ponent S. We normalize Fmm relative to the lowest observed 
Yukawa coupling, Fe, by defining Xmin = Fmin/Fe. The dark- 
est area is the Icr range, the second darkest area marks the 
2a range etc. 



where a product over all charged fermions i is under- 
stood. In comparison with Eq. ([2]), the gauge couplings 
have been dropped due to our focus on the primary con- 
straints due to the Yukawa couplings. The potential de- 
pendence on V in p(v,Ti...) is no longer present due to 
the assumption that the intrinsic probability distribution 
in V is roughly flat over our the allowed atomic window^. 
The atomic function ^(w, r„, F^, Fe) is summarized in 
Fig. [T] and p{Ti) is the probability distribution for the 
Yukawa couplings where we use (5 = 1, Tmin = 0.4Fe 
and Tfnax = 1.26. The normalization of L(v) is irrele- 
vant, we are only interested in estimating its shape. We 
calculate L(v) numerically by randomly populating the 
Yukawa couplings using the appropriate weight at differ- 
ent values of the vev. In particular we generate a set 
of three Yukawa couplings for the up-type quarks, for 
the down-type quarks and for the leptons. The smallest 
Yukawa coupling of each set is defined to be F„, Td and 
Fe respectively. The relative probability of satisfying the 
atomic constraints then yields the likelihood function. 

Let us briefly reiterate here the main assumptions 
which go into our analysis and explain the logic of how 
the application of a scale invariant probability distribu- 
tion for the Yukawa couplings can yield a constraint on 
the scale of the Higgs vev. First of all, we note that we 
extract the observed Higgs vev from measurements other 



L{v) = / dT, A{v,Tu,rd,Te) p(F,) (11) 



wc will return to these issues in our discussion of the uncertainties 
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FIG. 7: The likelihood function L for the Higgs vev as a 
function of v/vo (vo is the observed Higgs vev of 246 GeV), 
constructed with our favored scale invariant weight. The nor- 
malization gives the fraction of simulations which satisfied the 
anthropic constraints. 
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FIG. 8: This shows the same likelihood function L as in Fig. 
[71 but with the vertical axis on a linear scale. 



than the fermion masses, such as Mw- Now, our crucial 
assumption that the Higgs vev and the Yukawa couplings 
are statistically independent is used to infer the probabil- 
ity distribution for the Yukawa couplings independently 
of the Higgs vev. In practice, it essentially means that 
wc infer the scale invariant weight for the Yukawa cou- 
plings directly from the observed scale invariant weight 
for the quark and lepton masses for a fixed value of the 
Higgs vev. The probability distribution for the Yukawa 
couplings is taken to be a power law with an exponential 
close to the scale invariant case oi 5 = 1. For (5 = 1 it 



FIG. 9: The likehhood function L shown in FigQand Fig. [8] 
is shown here with both the horizontal and vertical a^ces on a 
linear scale. 



requires lower and upper endpoints. Whereas the lower 
endpoint of the distribution is inferred from the measured 
Yukawa couplings, for the upper endpoint the renormal- 
ization group quasi-fixed point has been used. Again, the 
assumption of statistical independence is crucial when we 
extract these endpoints and use them universally for any 
Higgs vev. Thus, we know that the Yukawa couplings 
have to not only be uniformly distributed on a log scale 
in a stretch extending over roughly six orders of magni- 
tude, but we know where on the log scale this stretch is 
located. Finally, a scale enters through the atomic con- 
straints on the light fermion masses, and together with 
the probability distribution of Yukawa couplings, it yields 
our constraints on the Higgs vev. Consider for example a 
large Higgs vev of 10^^ GeV. Since our Yukawa couplings 
from a scale invariant weight are required to lie roughly 
between 10"^ and 1, it would be impossible to get light 
fermion masses in the MeV range for such a large value 
of the Higgs vev. 

For our main result, we consider the case seen in Na- 
ture where only the m, d quarks and the electron fall 
within the anthropic window - all others quarks and lep- 
tons are heavier and should not be part of stable atoms. 
We discuss alternatives in the next section. We imple- 
ment this constraint by requiring that the second lightest 
up-type quark, the c quark, does not lie within the an- 
thropic window sketched in Fig. [T] when the rriu axis is 
replaced by rric, and analogously for the second lightest 
down-type quark. 

For our favored scale invariant weight, the result is 
shown as Fig. [71 [51 and [HI using log-log, log-linear and 
linear-linear coordinate axes. This is our primary result. 
We see that the distribution is peaked near the value 
vq observed in Nature and it extends over several orders 
of magnitude. The median value in this distribution is 
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FIG. 10: The likelihood function for different values of the 
lower endpoint in the Yukawa distribution, Vmin = OATe 
(solid) and Fmin = 0.04re (dashed). 

V — 2.25 Wo- The 2 — a range extends from 0.10 wo to 
11.7 vq. We observe that there is a steep upper cutoff in 
the allowed values of the vev which conies from the lower 
endpoint Tmin present in the scale invariant distribution 
of the Yukawa couplings. The relevant point where the 
likelihood function falls off faster than 1/t; is located at 
a few times vq. We conclude that vq would be a very 
typical value for the Higgs vev whereas values » 10 x wq 
would be very unlikely. 

7. UNCERTAINTIES 

In this section, we consider the changes in the likeli- 
hood function if we modify some of the features of our 
analysis. The two greatest effects conic from the varia- 
tion or removal of the lower endpoint in the fermion mass 
distribution, and from the possibility of extra quarks or 
leptons within the anthropic window. The first of these 
has the potential to significantly modify our results. 

First we can consider the changes if we use a differ- 
ent value of the lower cutoff. Changing the lower end- 
point by one order of magnitude within the context of the 
scale invariant weight produces the modification shown in 
Fig. [TUl The plot shows that for r™,„ = 0.04re the like- 
lihood distribution favors larger values of the Higgs vev. 
In this case, smaller allowed values of the Yukawa cou- 
plings are compensated by larger values of the Higgs vev 
to satisfy the atomic constraint. We see that the qualita- 
tive features of the spectrum remain. However, the value 
of V where the likelihood starts to falls off faster than 1/v 
is roughly 10 times higher if we divide Tmin by a factor 
of 10. We clearly need to address this issue of how much 
the most likely Higgs vev depend on the cutoff at low 
Yukawa couplings in the weight. 



10^ r- 

10'' - 

lO'i - 

10^ - 

1 - 



10-^ - 
10-' - 




lQ-6 I L , , , , , 1 

10-'' 10"' 10"^ 1 10^ 10' 10" 10** 
v/vo 

FIG. 11: The likelihood function resulting from a power law 
weight of exponent S — 0.86 without any lower endpoint Tmin- 

One major concern about the result of the previous 
section is that the shape of the likelihood function is 
determined at large values of the vev by the fact that 
the scale invariant weight has a lower endpoint to the 
Yukawa distribution. We address this issue by consider- 
ing a power weight with S less than unity, with no lower 
endpoint. Thus the power law behavior of the weight at 
low Yukawa couplings is not cut off at all and extends all 
the way down to zero. The Yukawa couplings can there- 
fore become arbitrarily small, without being in violent 
disagreement with the overall fermion mass distribution. 
Such a distribution has the potential to allow arbitrar- 
ily large values for the vev. In these cases, there would 
be situations where a high value of the vev is counter- 
balanced by very small Yukawas in order to satisfy the 
atomic constraints. 

In Fig. [Til we show the likelihood function that is ob- 
tained for S — 0.86 and Tmin — 0. While the maximum 
of the likelihood function remains close to the observed 
value we see that it never falls off faster than 1/v in the 
region studied which extends up to 10* x vq. That means 
the most likely values of the Higgs vev in this scenario 
would be very large. Because the power law weight is 
valid down to zero Yukawa couplings, we lose the con- 
straints for the Higgs vev to be in the neighborhood of 
the observed vq. While our analysis in Sec. [5] has shown 
that Tmin = is disfavored by over 2a this issue remains 
a serious caveat. It motivates further top-down studies 
of properties of the string theory landscape in order to 
identify the existence of a lower cutoff in the quark mass 
weight. 

Another uncertainty is the interesting possibility that 
more quarks beyond u, d, e fall in the atomic window. 
If we treat all the quarks independently as we have been 
doing we find that it is reasonably common that more 
quarks do fall in the allowed atomic window. For exam- 
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FIG. 12: The likelihood function allowing any number of light 
quarks in the anthropic window (dashed), along with that 
allowing u, d, e only (solid). 



pie, with the scale invariant weight the likelihood func- 
tion including extra quarks in this range is shown in 
Fig. [121 The allowance of extra small masses to fall in 
the atomic window makes the distribution peaked around 
smaller values in comparison to when we do not allow 
them to fall inside the window. The value of v where 
the likelihood starts falling off faster than 1 /v and there- 
fore the most likely Higgs vev remains almost unchanged. 
We do not see the disadvantage of this situation for the 
existence of atoms. 

However at this point it is useful to recognize a known 
flaw in our approximation that the Yukawa distributions 
of each flavor are treated as independent. Even if the 
original Yukawa couplings of the theory are distributed 
independently, the final output governing masses will not 
be independent. The original Yukawas exist in a 3 x 3 
complex matrix for each charge sector. The diagonaliza- 
tion of this matrix yields the final eigenvalues as well as 
rotation angles that go into the CKM and PMNS weak 
mixing matrices. It is well known that upon diagonal- 
ization, the eigenvalues of a matrix repel each other. In 
random matrix theory this leads to a repulsion of the final 
eigenvalues. For our case, this says that the u, d, e dis- 
tributions will be independent, since they come from dif- 
ferent charge sectors, hence different matrices. However, 
the likelihood of two quarks of the same charge falling in 
the allowed atomic range will be decreased by this repul- 
sion. We have studied this effect by generating random 
Yukawa couplings in a 3 x 3 matrix and diagonalizing. 
As discussed in in this case in order to approximate 
a scale invariant distribution of the final eigenvalues we 
start with an initial weight with S = 1.16, and we use 
Fmm — 0.4re. The resulting likelihood function is shown 
in Fig. I13[ where the solid curve only allows u, d, e in 
the anthropic window and where the dashed curve re- 
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FIG. 13: The hkelihood function resulting from a biunitary 
diagonalization of the Yukawa matrices allowing any number 
of light quarks in the anthropic window (dashed) , along with 
that allowing u, d, e only (solid). 



suits from allowing any number of light quarks in the an- 
thropic window. While the falloff at higher v is now less 
steep than in the result from diagonal simulations with- 
out matrix diagonalizations of Fig. [121 it is clearly falling 
off much faster than 1/v and the most likely value of v 
would be below 10 x vq. Comparing the result from the 
diagonalization of the Yukawa matrices in Fig. [131 with 
the corresponding result from diagonal simulations with- 
out matrix diagonalizations shown in Fig. I12[ we note 
that the two curves in Fig. [131 are closer together than 
the ones in Fig. [121 which is due to the repulsion of eigen- 
values in the matrix diagonalization case. 

Since our work is based on statistical independence of 
V and Ti, i.e. p{v,Ti) = p{v)p(Ti), there could be even 
more dramatic problems if there were correlations be- 
tween V and F;. For example, if in our weight p{Ti) both 
endpoints Tmin,Tmax were proportional to 1/v, there 
would be no anthropic constraints on v since varying v 
would keep the distribution of the masses the same. Be- 
cause we measure the Yukawa couplings at a single value 
of we cannot address this. 

We have mentioned previously the possibility that fur- 
ther anthropic constraints could come into play. These 
could change the detailed shape of the likelyhood func- 
tion because they could eliminate regions of parameter 
space which are somewhat different from our world. How- 
ever, since our parameters are clearly consistent with the 
constraints, the elimination of other regions would likely 
narrow the resulting likelihood function for v. This would 
change the shape of the function, but would not modify 
our basic conclusions about the compatibility of our value 
of V with the allowed anthropic range. 

One might also wonder if the likelihood function should 
also take into account the other great anthropic con- 
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straint, that on the chimping of matter in the universe 
to form stars and planets, which hmits in particular the 
cosmological constant. The likelihood function for the 
cosmological constant has been studied by Vilenkin and 
Garriga JJJ. The parameters of the Standard Model do 
of course also influence the cosmological constant. For 
example a shift in the value of the up quark mass (or the 
Higgs vev) by one part in 10'*° would shift the cosmologi- 
cal constant by 100%. However, for the cosmological con- 
stant to have an anthropic selection, the possible values 
of A must be densely packed and there should be other 
slightly different combinations of parameters that yield 
an anthropically allowed cosmological constant. For the 
rather narrow range of v that we probe, it seems very 
reasonable that the clumping of matter constraint has 
little influence on the likelihood of the Higgs vev. 

8. CONCLUSIONS 

The fundamental question that we are addressing is 
whether it is plausible that it is the atomic constraints 
that explain the low value of the Higgs vev in landscape 
theories. This is known to be the case at fixed values 
of the Yukawa couplings, but since these parameters also 
may vary in landscape theories, and since their magni- 
tudes seem to be peaked at low values, the answer is 
less obvious in a more general context. We have used 
experimental information on the distribution of masses 
to address this issue. We find that even if the Yukawa 
couplings are allowed to scan in a way that is favored 
phenomenologically, the likely values of the Higgs vev 
are close to the one observed in Nature. This supports 
the hypothesis that these constraints favor Higgs values 
similar to ours. 

This can be interpreted as a motivation for further ex- 
ploration of landscape theories. The value of the Higgs 
vev and that of the cosmological constant are two great 
"fine-tuning" problems of the fundamental interactions, 
and the presence of a landscape would change the way 
that we approach the issue of fine-tuning [15|. This is 
because both of these problems appear to have plausible 
resolutions through anthropic constraints that are appro- 
priate for landscape theories®. 

In [§[, it was argued that the anthropic constraints 
on the quark masses cannot be used to constrain the 
Higgs vev, by constructing a plausible scenario in which 
the weak interactions do not appear. In the context of 

* However, the "strong CP problem" is a fine-tuning puzzle that 
docs not appear to have an anthropic resolution |16|| . 



our exploration in this paper, this could be realized by 
taking w — > oo, F — 0, with their product fixed. This 
situation may arise in our scenario for power law weights 
with exponent S < 1 which extend down to zero Yukawa 
couplings. Such weights are disfavored by more than 2 — 
a in comparison to weights with a lower cutoff Tmin of 
the order of Fe as we infered from the measured ferniion 
spectrum. Nevertheless they are a serious caveat to an 
anthropic constraint for the Higgs vev. For the prefered 
weights with a lower cutoff Tmin of the order of Fe the 
most likely values are close to that seen in Nature. 

One of the strengths of our approach is that it does 
not rely on the ultraviolet completion of the fundamental 
underlying theory. Aside from our assumption of statis- 
tical independence, this input comes from the data on 
the quark and lepton masses. Both the issues of statis- 
tical independence and of weights without a lower cutoff 
could possibly be addressed in top-down studies of the 
landscape. 

The likelihood function that we have constructed is an 
estimator that tries to quantify the effect of the possible 
variation of the fundamental parameters on the range of 
allowed values for the Higgs vev. Further understanding 
of anthropic constraints may be able to narrow the like- 
lihood function further. Even if the range is narrowed, 
this is more of a consistency check than a prediction of 
landscape theories, since we already know that the ob- 
served value of the vev is anthropically allowed. However, 
it does provide further motivation for landscape theories 
and suggests that within these theories the hierarchy and 
fine-tuning problems associated with the vev are not as 
serious obstacles as they are in other theories. 
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